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We examine how to distinguish between unitary operators, when the exact form of the 
possible operators is not known. Instead we are supplied with "programs" in the form of 
unitary transforms, which can be used as references for identifying the unknown unitary 
transform. All unitary transforms should be used as few times as possible. This situation 
is analoguous to programmable state discrimination. One difference, however, is that the 
quantum state to which we apply the unitary transforms may be entangled, leading to a 
richer variety of possible strategies. By suitable selection of an input state and generalized 
measurement of the output state, both unambiguous and minimum-error discrimination can 
be achieved. Pairwise comparison of operators, comparing each transform to be identified 
with a program transform, is often a useful strategy. There are, however, situations in which 
more complicated strategies perform better. This is the case especially when the number of 
allowed applications of program operations is different from the number of the transforms to 
be identified. 

Keywords: unambiguous discrimination; optimum discrimination; operator comparison; 
generalized measurements 



1. Introduction 

The problem of discriminating between operators in quantum mechanics is closely 
related to that of discriminating between quantum states dj). In order to distin- 
guish between operators, a reference state is transformed by the different operators, 
and then measurements are performed on the result, so that one is ultimately dis- 
tinguishing between the possible output states. Discriminating between two known 
operators has already been treated pj-Q). Reconstructing a single unknown oper- 
ator from measurements on a series of states that have been transformed by it is 
known as operator, channel, or quantum process tomography Q). A final problem, 
operator comparison, which determines whether two unknown unitary operators 
are the same or different (0, H), is closely related to the problems considered here. 

Recently state discrimination problems have been considered in which not all 
of the states are known (fi^ -[73). One may know some of the states, but is only 
provided with examples of the others. For example, one is given a particle that is 
guaranteed to be in one of the states \ip\) or | -i/>2 ) ■ The states \ipi) and |^) are 
unknown, but instead, an example of each of these states is provided. The object 
is to use the reference states to best determine whether the given particle in the 
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unknown state is either or IV^)- Another scenario is that |V>i) is known, but 
|^2) is represented by a reference particle, and one is presented with a particle 
guaranteed to be in either or | -i/>2 ) • Again the task is to determine the state of 
the latter unknown particle. 

In this paper we address problems in the realm of operator discrimination, which 
are the natural analogues of those, decribed above, for quantum state discrimina- 
tion. Instead of having a physical description of the operations, they are given as 
actual operations (boxes) which we are allowed to apply to input states. The situa- 
tion is reminiscent of oracle quantum computation (|Xj) , in which a key component 
is a device that carries out one of a number of possible transformation. Generally, 
we would like to decide which particular unlabelled boxes match/correspond to 
labelled reference boxes or in which order the unlabelled boxes are placed. The 
problem can be viewed as programmable pattern matching, "programmable" be- 
cause we do not know the actual operations as they are only supplied as reference 
"programs" . 

We shall consider both minimum error, and unambiguous discrimination tasks. 
In minimum-error discrimination, the object is to always make a guess and to 
minimize the probability of making a mistake. In unambiguous discrimination, the 
procedure is allowed to fail sometimes (that is, we are allowed sometimes to refuse 
to make a guess) but we never want to make a mistake if we do decide to make a 
determination. The object is to minimize the probability of failure. 

We wish to examine whether good procedures for these decision problems can 
be constructed from pairwise operator comparison. The task is to compare two 
unknown operators in order to determine whether or not they are the same. This 
problem was considered in Q), and the best procedure found there was to use a 
singlet state (Q). Each operator is applied to one of the two qubits in the singlet and 
the result is then measured to see if it has a component in the symmetric subspace, 
upon which the operators must have been different. We would like to see if this 
procedure can be used as a basic unit for more complicated operator comparisons 
as discussed above. We shall see that there are cases where singlet-based pairwise 
comparison strategies are not optimal. 



2. Pattern matching for unitary operators 

The general scenario which we consider is as follows: Alice/Bob are given reference 
boxes, labelled either U or V, which implement unitary operations on single qubits. 
The forms of these two unitary operators are not known to any of the parties. 
Charlie is given a number of boxes which are unlabelled, but each is guaranteed 
to perform either U or V. Our task is to determine what boxes Charlie possesses 
by applying these boxes to a multi-qubit reference state and then measuring the 
results. We shall assume that each box can only be applied once and that U and V 
are equally and independently likely to be any qubit unitary. A variety of scenarios 
is possible, depending on which boxes are available for Alice, Bob and Charlie. 

2.1. One reference box and one unlabelled box 

We begin with the simplest case of one reference box and one unlabelled box which 
may of may not match the reference (Fig. [T^)- Alice has a box that performs an 
unknown operation U and Charlie has one that eitherperforms U or an unknown 
operator V. This is similar to operator comparison 0) where it was found that 
making use of the two-qubit singlet state, \i/j~)ab = (|0)a|1)b — 10)^1 1}b)/v2, by 
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Figure 1. One unlabelled box. a) Alice has a reference box U, and Charlie has an unlabelled box which 
either performs U or a different unknown unitary V. b) Alice and Bob each have a different reference box 
implementing U and V respectively, and Charlie either has a box which performs U or V. We prepare a 
test state |<ip) which is acted upon by the operators and is finally measured leading to an outcome IIj . 

sending one qubit through each box is a good strategy. After the qubits have been 
acted upon by the boxes, we measure the result to see if it lies in the symmetric 
or antisymmetric subspace. If the state is found to be symmetric, then the two 
unitary transforms must have been different. 

This is in fact the optimal strategy also for the situation we are now considering, 
if we confine ourselves to two qubits and demand an unambiguous result. We have 
to apply both transforms to a test state |^), giving two possible cases 



\^i)ac = U® U\^)ac and \^ 2 }ac = U <g> V\& 



'AC- 



(1) 



The forms of the single-qubit operators, U and V, are unknown, and so we must 
average over U and V to obtain the possible output states 



Pi = J du\^x) AC (^i\ and p 2 = j dudv\$> 2 ) Ac \ 



(2) 



between which we should distinguish. The integrals are Haar integrals over SU(2), 
and can be performed using the following parameterization. We first set U = 
exp(— i#e • a/2), where e is a unit vector equal to (sin \i cos (p, sin \i sin (ft, cos /jl), 
with < 6 < 2tt, < (f> < 2ir, and < \i < ir, and a = {a x , a y a z } is a vector of the 
usual Pauli matrices. The second operator, V, is parametrized similarly. We then 
have that 



du 



1 

4^ 



sin 



sin fj,d6d(f>dfj,. 



(3) 



A general 2-qubit test state can be written 

\V)ac = M~) + b\^ + ) + <#~) + 



(4) 



where 1^) = (|01) ± \10))/V2 and 1^) = (|00) ± |11))/V2 are the usual Bell 
states, and \a\ 2 + \b\ 2 + |c| 2 + \d\ 2 = 1. Performing the integrals for this initial test 
state, see Appendix [A] we find that 



p x = \ a \ 2 P a + \(\b\ 2 + \c\ 2 + \d\ 2 )P s 
92 = jIa ® Ic- 



(5) 
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Here P a and P s are the projectors onto the antisymmetric and symmetric subspaces 
respectively. As the support of p2 is all of the two-qubit space, p\ can never be 
unambiguously identified, as its support will always be contained in that of p%. To 
unambiguously identify the state p2, the initial test state should lie entirely either 
in the antisymmetric or in the symmetric subspace, and we look for a component 
of the final state in the symmetric or antisymmetric subspace, respectively. 

It is straightforward to check that starting with a singlet state by choosing \a\ = 1 
will give a higher success probability than starting in the symmetric subspace. The 
positive operator valued measure (POVM) elements for the optimal measurement 
are the projections Ilo = P a , corresponding to the inconclusive result (failure), and 
II2 = P s , corresponding to identifying the state as pi. Therefore, assuming that p\ 
and P2 are equally probable, the probability of the measurement failing is 

p f = ^[Tr(riopi) + Tr(n p 2 )] = |, (6) 

so that the success probability is p s = 3/8. An initial test state in the symmetric 
subspace would give a success probability of p s = 1/8. 

Now let us consider the minimum-error strategy. The optimal minimum-error 
measurement for distinguishing between two states p\ and P2, with prior probabil- 
ities p\ and P2 , has the error probability (fTj ) 

1 

Pe = ~ 2 



1 - TV (PiPi -P2P2) 2 , (7) 



where we shall consider the case p\ = P2 = 1 /2. Substituting in the density matrices 
from Eq. ([5]), we find that 



TiV(pi - P2 f 



1 I |2 



(8) 



This is clearly maximized, and the error probability minimized, when \a\ 2 = 1, 
which occurs when the input state is a singlet. In that case we have that IIi = P a 
corresponds to Bob's box being U, and H2 = P s corresponds to Bob's box being 
V, and the probability of error is 1/8. For the singlet test state, this is the optimal 
measurement, as the optimal minimum-error measurement. 

Note that the unambiguous discrimination strategy only gives us an answer when 
Charlie's unknown box is V. If Alice instead has a reference box V, we can get an 
answer only when Bob's box is U. The measurement and success probability stays 
the same. 



2.2. Two reference boxes and one unlabelled box 

Let us now consider a decision problem where Alice has a reference box U, Bob 
has reference box V, and Charlie has an unlabelled box guaranteed to be either U 
or V (Fig. [lb). First let us look at a three-qubit test state \*$>)abc and apply the 
boxes held by Alice, Bob and Charlie to their respective qubits. A general pure 
three-qubit state can be written as 



^)abc = ao|000) + oi|001) + . . . + a 7 |lll), 



0) 
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with Y, N 2 = 1. This means that we are trying to discriminate between the states 

\*i)abc = U ®V ®U\^) ABC (10) 
\^2)ABC = U® V® V\^)abc- (11) 

Averaging over U and V means that we have to distinguish between the states 

Pi = J dudvU ®V ® U\^)abc{^\U ] (g> V ] <g> U ] 

P2 = J dudvU ®y® V\^)abc(^\U^ <g> (8) (12) 

Again using the results in Appendix A, we get 

Pi = \p%c ® ^b(|oi - a 4 | 2 + |a 3 - a 6 | 2 ) 



+^lc "X 1 ^b(|oo| 2 + Kl 2 + |ai + a 4 | 2 + \a 3 + a 6 | 2 + |a 5 | 2 + |a 7 | 2 ) 
P2 = \pbc ® /a(|oi - a 2 | 2 + |a 5 - a 6 | 2 ) (13) 
+\ p bc ® /a(|oo| 2 + |o 3 | 2 + |ai + a 2 | 2 + |a 5 + a 6 | 2 + |a 4 | 2 + |a 7 | 2 ). 



For unambiguous discrimination to be possible, the support of either p\ or p 2 or 
both has to be less than all of the three-qubit space. Without loss of generality, 
let us assume that p\ has support in either P%c ® Ib or P^ c (g) Ig, but not in 
both. The first possibility clearly corresponds to a test state which is a singlet in 
qubits A and C, tensored with any state for qubit B. This strategy has an overall 
success probability of 3/8. The optimal measurement tests whether or not we have 
a singlet in qubits A and C, and if not, this unambiguously identifies the state as 
pi- We could also start with a singlet state in qubits B and C, and test whether 
we still have a singlet in these qubits after application of the unitary transforms. 

The other possibility is that p\ has support in the space symmetric in qubits A 
and C. In this case p2 either has support in all of Hilbert space, or only in the 
space symmetric in qubits B and C. The latter alternative gives a higher success 
probability, as we then have a nonzero probability to unambiguously identify both 
pi and p2- We then have 

Pi = \p S ac ® Ib and p 2 = ^P S BC ® I a- (14) 

The POVM that unambiguously discriminates between these two density matrices 
must have the element that detects p\ be proportional to P^ c ®Ia, because the only 
vectors orthogonal to p 2 are of the form \i[>~) bc ® \v)a, where \t])a is an arbitrary 
state in the space of qubit A, and by symmetry, this POVM element should treat 
all vectors in the space of qubit A equally. Similar considerations apply to the 
element that detects p 2 . The only remaining task is to determine the constant of 
proportionality that is allowed by the requirement that the third POVM element, 
which correspond to the failure of the measurement, be positive. Doing so, we find 
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that the optimal measurement in this case has measurement operators 

2 



III 

n 2 
n ? 



3 Pbc^Ia 

\PAC®lB 



I -III 



n 5 



(15) 



and the overall success probability is 1/6. This is higher than 1/8, which is the 
result if we choose to randomly test for a singlet either in qubits AC or qubits BC. 
But it is still much lower than 3/8, which we achieve if we start with a singlet state 
either in qubits AC or qubits BC. 

For the minimum-error task, we can also try to use a state that approximates a 
singlet state in qubits AC and BC as well as possible, without obeying the strict 
symmetry rules necessary for the unambiguous discrimination strategy. The state 

\^)abc = ^(l 10 °) + I 010 ) + I 011 ) + 1101)) - ^=(1001) + |H0)) 

has the property that its reduced density operators pac an d Pbc have overlap with 
the singlet state 

(?P~\PAc\ip~} = (iP~\pbc\iP~} = 3/4 



where \ip~) is a singlet. It is straightforward to prove that the fidelity of 3/4 is the 
best possible if one demands that the fidelities in qubit pairs AC and BC should 
be equal. After Alice, Bob and Charlie have applied their boxes to obtain (|10p and 
(fTTT) . and then averaging over U and V, the two possible output states are 



Pi 



3 1 

pa i ids 

S AC + 24 AC 



Ib, 



P2 



P%C + 



1 

24 



BC 



Ia, 



which then have to be distinguished. The optimum minimum error (Helstrom) 
strategy for distinguishing between two density operators p\ and p2 has the error 
probability 



Pe= 2< 1 



(Pi/01 -P2P2Y 



1/2 



1 



where Aj are the eigenvalues of pipi — P2P2- In our case pi = P2 = 1/2, and 
1 " 



PiPi - P2P2 



3 1 

"pa , ps 

V AC + 24 AC 



'•J - ( lp a BC + 



Ia 



\Plc + Y lAC ) ® Ib ~ \l P BC + Y lBC ! ; ' 1 



-[P%C®lB-P a BC ®lA] 

(|100)(100| + |011)(011| - |010)(010| - |101)(101| 
— 1100> <001| - |110)(011| + |010)(001| + |110)(101| + h.c.) . 
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This operator is six-dimensional, but happens to be block diagonal with two blocks 
of size 3x3. This is easily seen since the off-diagonal elements only couple states 
with an equal number of 0's and l's. Ordering the basis vectors within the blocks 
{|010),|100>,|001>} and {|101>, |011), |110)} , each of the blocks takes the form 




The blocks have eigenvalues ±l/(4\/3) and zero, and the sum of the absolute 
values of all eigenvalues is therefore 1 / \/3. This gives an error probability of P e = 
(1 — l/\/3)/2 ~ 0.211325 which is higher than 1/8. Therefore, the simple singlet 
strategy (a singlet either in pair AB or AC) does better. 

2.2.1. Doing better than pairwise operator comparison 

So far, the best strategies we have found are based on pairwise operator compar- 
ison. Can we do better? It turns out that if we use a four-qubit state then we can. 
In particular, we shall show that it is possible to improve the success probabilities 
for unambiguous discrimination. We first note that we can get distinct signals for 
both operators in the case of unambiguous discrimination if we use four qubits. Let 

Wabcd = ^(|YOac|0)b|0) d + \i/>-)bc\0)a\1)d). (16) 

As usual, Alice, Bob and Charlie apply their respective boxes but no operation is 
applied to the fourth qubit D. It is best to express the resulting states as density 
matrices. If Charlie's box performs U, the operator U (&V ®U ® I is applied to 
\^}abcd, and after averaging over both U and V we obtain 

Pi = \Plc ®Ib® (\0)d + \\1)d) ( d(0| + \ D (l 

+^P%c®Ib®\1)d(1\. (17) 

If Charlie's box performs V, the operator U®V®V®Iis applied to \^)abcd, 
and after averaging over both U and V we obtain 

P2 = \P%C ®Ia® (\1)d + ^|0)d) (d(1| + \d(0\ 

+^Pbc®Ia®\0)d(0\. (18) 

We will consider only the case of unambiguous discrimination and first choose the 
POVM elements to be 

111 = P s bc® !a® 

11 2 = px c ®i b ®\o) d (o\ 

n = i-n 1 -n 2 = PZ c ®iA® |i)d(i| + p% c ®Ib® |o) d (o|. (19) 



The element corresponds to Charlie's box being U, II2 corresponds to it be- 
ing V, and Hq corresponds to the indeterminate result. We find that Tr(IIip2) = 



September 23, 2011 0:12 Journal of Modern Optics operatordecisionJModOptArxiv 



8 Andersson, Barnett, Hillery and Oi 




Figure 2. Two references and two unlabelled boxes, a) b) c). 



Tr(Il2/0i) = 0, which guarantees that the discrimination is unambiguous. In addi- 
tion, we have that Tr(IIipi) = Tr(Il2P2) = 3/8, which, assuming that the boxes for 
U and V are equally probable, gives a success probability of 3/8. This is the same 
success probability as with the simpler two-qubit procedure, but we now have a 
probability of unambiguously identifying both boxes, and not just one. 

So far, in terms of the overall success probability, we have not gained anything. 
We could in fact also have achieved the same result using a two-qubit strategy, 
by choosing to compare Charlie's box against either Alice's U or Bob's V with 
probability 1 jl each. The fourth qubit above acts like a quantum coin. By using the 
four-qubit test state in (fl~6]) we can, however, do better than the success probability 
of 3/8. Though this test state is most probably not optimal, we can try to improve 
the POVM. We make use of the results of subspace discrimination to find the 
optimal POVM for the discrimination of the ranges of the two possible density 
matrices While this does not give the optimal POVM for discriminating 

the two density matrices, it does give one that is better than (TT9l) based on two- 
operator comparison. The details of the POVM that optimally discriminates the 
ranges of p\ and p2 are given in Appendix B. One finds that each density matrix 
can be detected unambiguously with a probability of ~ 0.43, which is greater than 
3/8 = 0.375. Therefore, we can conclude that there are better procedures than 
simple pairwise operator comparison for solving our three-box decision problem. 
The optimal strategy remains to be determined. 



3. Two reference boxes and two unlabelled boxes 

When Charlie possesses two unlabelled boxes, there are several scenarios we can 
consider (Fig. [2]). Alice and Bob can possess U and V respectively, Charlie is given 
two boxes either both U or both V, and the task is to determine which boxes 
Charlie has. We can also consider the case in which Alice and Bob possess U and 
V respectively and Charlie is given two unlabelled boxes, say C\ and C2, guaranteed 
that either C\ = U and C2 = V or vice versa, and our task is to determine the 
order in which U and V appear. For this task, we will also consider the case where 
the reference boxes that Alice and Bob are given are either both U or both V. 



3.1. Two reference boxes and two of the same unlabelled box 

Let us consider the case of Alice having box U, Bob box V, and Charlie a pair of 
identical boxes, both of which are U or V (Fig. [2^). A good strategy here is to use 
four qubits in the form of two singlet states, \^)abc 1 c 2 = \'4>~)ac 1 \'4 , ~)bc 2 an d 
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apply the boxes. If both boxes C\ and C 2 are V, then the output state is 

\*i)ab Ci c 3 = {U®V)\ii-)A Cl (V®V)W-) B c 2 , (20) 

and if both C\ and C 2 are U, then we have 

\*2}ab Ci c 2 = (U® U)\^-) ACl (V U)\^-) BC2 . (21) 

If we now average over both U and V, we obtain the possible output density 
operators 

Pi = \lA Cl ® Pbc 2 

P2 = \p%c,®Ibc 2 . (22) 

between which we must distinguish. 
Let us first consider unambiguous discrimination. The best POVM elements are 

n 1 = PX Cl ®P£ C2 
n 2 = p a ACi ®p s BC2 

no = p s ACi p s BC2 + P a ACl ® P a BC 2 , (23) 

where IIo is associated with the failure outcome. We find that, as is necessary, 

Tr(nip 2 ) = Tr(n 2 pi) = 0, and that Tr(U lPl ) = Tr(U 2 p 2 ) = 3/4. With both 

situations being equally likely, our overall success probability is then 3/4. 

For minimum-error discrimination, we always need a result; the measurement is 

not allowed to fail. This means that we have to do something with the operators 

appearing in IIo. We first note that the outcome corresponding to P% c <8> P%c 2 

never occurs, so it can be added to either IT or 1T2 without any effect. The outcome 

corresponding to P%q ®P b 2 can occur, and for either density matrix does so with 

a probability 1/4. When it does, we have no idea which density matrix we have, 

so we have to guess. One way of guessing is that we can always guess p 2 , and we 

(2) 

shall be wrong half the time. This corresponds to adding P%q x ® P]3C 2 ^° (^ or 
equal prior probabilities of p\ and p 2 , it does not matter how we guess, but if the 
probabilities are unequal, we should always guess the more likely one.) Therefore, 
for minimum-error discrimination we can choose for our POVM operators 

n 2 = p% Ci ® p BC2 + p% Ci ® p% C2 , (24) 

and the probability of error will be 1/8. 



3.2. Two reference boxes and two different unlabelled boxes 

We now turn our attention to a variant of our decision problem in which the object 
is to determine which box corresponds to which operator. Charlie has two boxes 
C\ and C 2 , and either box C\ performs U and box C 2 performs V, or vice versa. 
Our object is to determine which of these two situations holds. This version of the 
problem is sometimes equivalent to the version with one unknown box considered 
in previous sections. For example, if in addition to Charlie's boxes, there is only 



September 23, 2011 0:12 Journal of Modern Optics operatordecisionJModOptArxiv 



10 Andersson, Barnett, Hillery and Oi 

one program box U, we now have to distinguish between U® V ® U and V ® U® U . 
As we are averaging over U and V, we can however interchange these in the second 
density matrix, so that V ® U ® U becomes U ® V ® V. This scenario is equivalent, 
therefore, to the third box being unlabelled, with one reference box each of U and 
V. All the strategies we have discussed above that apply to this three-box situation 
can therefore be carried over to the order-finding version. It remains, therefore, to 
consider the four-box situations, in which Alice and Bob both have reference boxes.. 

3.2.1. Two reference boxes of U 

In the case in which Alice and Bob both have the same reference box, say U, and 
Charlie has either C\ = U and C2 = V, or vice versa ((Fig. O))), we can map this 
case directly onto the case in Section 13.11 We note that the two possibilities for 
ABC1C2 are either U ®U ®U ®V or U ®U ®V ®U . As we have to average over all 
possible U and V, this is equivalent to trying to distinguish between U ®V ®U ®U 
and U ® V ® V ® V where we have simply re-ordered the boxes (C1C2AB) and 
interchanged U and V in the second case. 

3. 2. 2. One reference box of U and one reference box of V 

One might wonder if replacing one of the U operators by a V in this scheme 
would improve the results. That is, suppose that Bob has a box that performs V 
instead of U (Fig. [2b). If box C\ is U and box C2 is V, then we now have the state 

\Vi)ab Ci c 2 = {U®U)\i>-) ACl {V®V)\i>-) BC2 , (25) 
and if box C\ is V and box C2 is U, we have 

\^>2)abc iC2 = (U®V)\ip-) ACl (V®U)\ip-) BC2 . (26) 
We average these two states over U and V giving us the two density matrices 

Pi = |*l)ABC 1 C 2 (*l| 

92 = \{\p s b Ci ® Pac 2 + P a B Cl ® Pac 2 )- (27) 
In comparing these two density matrices, it is useful to express \^>)ABdC 2 as 

\*)AB Cl C 2 = \{\^ + )BcM + )AC 2 + \^-)BcM-)AC 2 -mB Cl AC 2 -\H) B Cl |00) AC 2 )- 

(28) 

The first thing to notice is that the support of p\ is contained in the support of 
P2, so that if we are considering unambiguous discrimination, we will be able to 
positively identify p2, but we will never be able to positively identify p\. 

We choose II2, the POVM element corresponding to P2, to be the sum of three 
projection operators. The first is the projection onto the part of 'H s BCl ®Hac 2 
is orthogonal to the vector 

\q) = H + )bcM + )ac 2 ~ \00) BCl \n) A c 2 - \ll) BCl \00) A c 2 ), (29) 

the second is the projection onto rl s BCi ®'rV A o 2 , and the third is the projection onto 
rl a BCi ® T~i s A c 2 ■ The notation here is that rl s - k is the symmetric subspace of qubits 
j and k, while Hj k is the antisymmetric subspace of qubits j and k. The POVM 
element corresponding to failure is IIo = I — II2. We find that T^IT^/^) = 2/3, 
which means the our overall success rate is 1/3; we succeed with a probability 
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of 2/3 if we have p2, which happens with a probability of 1/2, and we get an 
indeterminate answer if we are given p\. 

We have found that replacing one of the U operators by V has made the situation 
worse. This makes sense as detecting a difference in either the pair AC\ or BC2 will 
identify the order of U and V for boxes C\ and Ci- If we match up identical oper- 
ations in both singlet pairs, then we have no chance of unambiguously identifying 
the order of boxes C\ and Ci- Matching different transforms up in both pairs gives 
a higher probability to detect a difference, as opposed to when there is a difference 
in only one singlet pair. This cannot, however, make up for our inability to detect 
a difference when identical transforms were matched up. It may be that in this 
situation, a strategy with more than four qubits would help. This is because the 
situation is somewhat similar to the three-box situation in Section \2. 2. li There, we 
would like to compare the unlabelled box against a reference box which does not 
match it, in order to unambiguously identify the unlabelled box. We cannot know 
which reference box to choose, however, and half the time we will choose to com- 
pare the unlabelled box against an identical reference box. The four-qubit strategy 
somehow helps us beat the odds. Here, we would like to compare both unlabelled 
boxes with reference boxes which do not match, but half the time we make the 
wrong choice and match the unlabelled boxes up with matching reference boxes. 
A strategy with more qubits might do better, but this is still an open problem. 



4. Conclusion 

We have discussed a number of operator decision problems. In each we are set the 
task of determining whether boxes perform the same or different unitary transfor- 
mations, U or V. We have no prior information of the forms of the transformations, 
but, instead, we have reference boxes that are known to perform the transforma- 
tions U or V. 

We began by considering one reference box and one unknown box. This situ- 
ation reduced to operator comparison, and for both unambiguous discrimination 
and minimum-error discrimination it is best to send into each box one of the par- 
ticles from a singlet state. We then went on to examine the situation where one 
has two reference boxes and one unknown box. This can be thought of as a kind 
of pattern matching. Each of the reference operators represents a pattern, and we 
are trying to determine which of the two patterns the unknown operator matches. 
We studied both unambiguous and minimum-error strategies based on pairwise 
operator comparison. We also studied a minimum-error strategy based on an en- 
tangled three-qubit state that best approximates two singlet states. This strategy 
was not as good as pairwise operator comparison. Using subspace discrimination, 
we found an unambiguous strategy that is better than pairwise operator compar- 
ison, which demonstrates that the pairwise comparison strategy is not optimal. 
Finally, we studied situations with two reference boxes and two unknown boxes. 
For these cases, we only studied strategies based on pairwise operator comparison. 

Our objective was not to present a comprehensive theory of operator compari- 
son problems, but to to gain some idea of the issues involved and to present some 
strategies that will work for specific problems. One conclusion is that pairwise op- 
erator comparison works well, but is not always optimal. Finding optimal strategies 
is a subject that deserves further study. 
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Appendix A. 

Here we present some additional information on the Haar integrals that were used 
extensively throughout the paper. As was noted, for SU (2) these integrals can be 
performed using the following parameterization. We first set U = exp(—i0e- <?/2), 
where e is a unit vector equal to (sin/xcos<^>, sin/xsin<^>, cos/x), with < < 2ir, 
< 4> < 2-7T, and < fx < n, and a = {cr x ,a y ,a z } is a vector of the usual Pauli 
matrices. We then have that 

du = -—r sin 2 ( - ] sin fid9d(j)dfi. (Al) 



4vr 2 V 2 . 
By explicit integration we find that 

J dui/|0><0|l7t = J du(J\l){l\rf = \j, 

J duU\0){l\U ] = J du(J\l){0\U* = 0. (A2) 

The unitary invariance of the measure then implies that for any qubit state 

j duU\Tl>)ty\rf = \h 

J duU\ip){ip±\U^ = 0, (A3) 

where (ip\^±) = 0. 

Now let us move on to two qubits. One can explicitly evaluate the integrals 
J duU ® [7|00)(00|C/ t ®[/t = ^P s 
J duU ® C/|00)(ll|C/ t = 

J duU ® U\i/)-){00\rf ® £/ f = 0. (A4) 
The unitary invariance of the measure allows us to immediately extend these results 

J duU ®U\^ S ){^ S \U^ ®t/ f = ^P S 
J duU ® U\ip s ){ip sl _\U^ <g> C/ f = 

J duU ® U\i)-){tp s \U ] <g> U ] = 0, (A5) 
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where \ip s ) and \ip s ±) are in the symmetric subspace, (Tp s \tps±) = 0, and \ifi a ) is in 
the antisymmetric subspace. 



Appendix B. 

Here we present the details of the POVM that optimally discriminates between 
the ranges of p\ and p 2 in Eqs. (fT7|) and (fT8|) . Let S\ denote the range of p\. It is 
spanned by the orthonormal vectors 

= -^=Hj-)ac\0)b (\0)d + i|l>uj , |«b) = |H)ac|0)b|1)d 
\u 2 ) = A|^-) AC |i) fl ^| ) D + l\l) D J , \ U6 ) = |11) AC |1) B |1) D 

|«3> = |00)ac|0)b|1)d, \u 7 ) = |V + )ac|0)b|1)d 

|«4> = |00)ac|1)b|1)d, \u 8 ) = \^ + ) ac \1)b\1)d- (Bl) 

Here, \if) + ) = (|0)|1) + |l)|0))/\/2- In this basis, we have 

5 1 8 

pi = ^d u i)^il + KX^I) + Ye 5^ l^^'l- (B2) 

Now, let S2 denote the range of p 2 . It is spanned by the orthonormal vectors 
h) = -^|^-)bo|0)a + \\V)Dj , |«5> = |H)bc|0)a|0)d 

\v 2 ) = ^=\^-) bc \1)a (\1)d + ^\0)d^ , k) = |H)bc|1)a|0)d 

|«3> = \00)bc\0)a\0)d, \vj) = \^ + )bc\0)a\0)d 

|«4> = |00) bo |1)a|0)d, K> = \^ + )bc\1)a\0)d- (B3) 

In this basis 

5 1 8 

P2 = t^(\vi)(vi\ + \v 2 ){v 2 \) + T^J2 \ v 3K v i\- ( B4 ) 

i=3 

Examining these vectors, we see that \u^) and \uq) are orthogonal to S 2 , and 
|z;3) and 1 1*4) are orthogonal to S±. The remaining vectors can be arranged into 
two sets of overlapping subspaces. The subspaces spanned by {|wi), |it4), \u?)} and 
by { \v± ), ^4), \v?}} overlap and the subspaces spanned by {\u 2 ) , \u^) , \us)} and 
{^2)) 1^5), 1^8)} overlap. The subspaces in the first set are orthogonal to those 
in the second. 

We can now apply the results of reference (T^ ) to discriminate the two three- 
dimensional subspaces and to discriminate the two two-dimensional subspaces. 
These results are then combined with those of the previous paragraph to give 
us the two optimal POVM's for unambiguously discrimination Si and S 2 . In order 
to exhibit these POVM's it is necessary to define a number of additional vectors. 
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We have in S± 



\ u 2) 



-L(|U4) + V2|U7» 

-^(|u 5 ) + V2|u8)) (B5) 



and in S2 



\ v 2) 



A=(\ V4 ) + V2\v 7 )) 

-L(\v 5 ) + V2\v 8 )). (B6) 



Both \ui) and 1 7x2) are constructed so that they are orthogonal to S2 and |%) and 
|i)2} are orthogonal to S±. Our discrimination problem is now reduced to discrimi- 
nating between two sets of two-dimensional subspaces. The first set consists of the 
subspaces spanned by {|i*i), 1%)} and 1^3}}, and the second set consists of 

the subspaces spanned by {| 1*2)1 1 ^4)} and {\v2}, $4)}, where 



\ v 3) 



^(V^|«4> " |«7» |U4) = ^(^K) " |«8» 

= -^(V2|«4> - \vr)) \v 4 ) = A=(V2\vb) - \v s )). (B7) 



The two sets are orthogonal to each other. The next step is to find the Jordan bases 
for each set of subspaces (fT^h . In each case, we find that the two subspaces have 
a vector in common. Throwing these vectors out, because the clearly do not help 
discriminate the subspaces, we are left with two sets of one dimensional subspaces 
to discriminate; \u^) and \vs) on the one hand, and \uq) and \vq) on the other. 
These vectors are defined by 



\ v 5) 



^=(Vs\ui - V5|«3» |«6) = ^^(v3|u2 + V5|iU)) 

- VE\v 3 )) \v 6 ) = tA^K + VE\v A )). (B8) 



Our final step is to define the vectors 



5 3 5 3 

\V5) = t(I^5) + 7^5)), \zb) = 7(1%) + 7^5)) 



4 m "/ 5 ' 

|2/6) = 7 = (I*%) + 7 1 K)) l«6> = 4(l«6>+gl«6))- (B9) 



These vectors have the property that (yj\uj) = and (zj\vj) = 0. for j = 5,6. 
The POVM operators for optimally unambiguously distinguishing Si and S2, and 
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thereby distinguishing p\ and p2, are 

III = \u 3 ){u3\ + \uq)(u 6 \ + \ui){ui\ + \u 2 ){u 2 \ 
5 

+ g(l«5)(«5| + 

n 2 = |W 3 )(U3| + |^6> {^6 1 + \vi){vi\ + \V2){V2\ 

+I(b 5 )(y 5 | + |y6)(y6|) (Bio) 
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